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Abstract 

We present a semi-perturbative approach which yields an infrared-stable fixed point in the 
Ginzburg-Landau for N = 2, where N/2 is the number of complex components. The calculations 
are done in d = 3 dimensions and below T c , where the renormalization group functions can be 
expressed directly as functions of the Ginzburg parameter n which is the ratio between the two 
fundamental scales of the problem, the penetration depth A and the correlation length £. We find a 
charged fixed point for k > l/y/2, that is, in the type II regime, where Ak = k — l/y/2 is shown to 
be a natural expansion parameter. This parameter controls a momentum space instability in the 
two-point correlation function of the order field. This instability appears at a nonzero wave- vector 
Po whose magnitude scales like ~ Ak^, with a critical exponent (3 = 1/2 in the one-loop approxi- 
mation, a behavior known from magnetic systems with a Lifshitz point in the phase diagram. This 
momentum space instability is argued to be the origin of the negative 77-exponent of the order field. 

PACS numbers: 74.20.-z, 05.10Cc, 11.25.Hf 
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I. INTRODUCTION 



Whereas the field-theoretic understanding of the universal critical properties of the su- 
perfluid phase transition has become quite satisfactory |l], |2|, the renormalization group 
(RG) have so far not been completely successful for the superconducting phase transition. 
The standard field theory investigated in this context is the Ginzburg-Landau (GL) model, 
initially in 4 — e dimensions by Halperin, Lubensky, and Ma (HLM)|3[], who generalized a 
calculation of Coleman and Weinberg |Q in four dimensions, and found no infrared stable 
fixed point at the one-loop level. The renormalization group (RG) flow was shown to run off 
to infinity in the g — /-plane if N < 365.9, with N/2 being the number of complex compo- 
nents of the scalar order parameter field. Here g and / are the dimensionless versions of the 
quartic self-coupling u of the order parameter and the square charge e 2 , respectively. HLM 
interpreted this result as indicating that the phase transition in a superconductor - which 
corresponds to a order parameter with N = 2 - is of first-order. They supported this scenario 
by a mean-field estimate of the effective potential obtained by integrating out the vector 
potential, which led to a term || ~ |?/'| 4_e , which leads to first-order phase transition. This 
was in contrast to Monte Carlo simulations of lattice superconductors in 1981 by Dasgupta 
and Halperin || where a clear second-order phase transition was found. It also contradicted 
work by Lawrie and Athorne |7J who performed a RG study in a non-linear realization of the 
GL model and found no indication of a first-order phase transition. Moreover, the first-order 
scenario was also contradicted by experiments on the smectic-A to nematic phase transition 
in liquid crystals || which according to de Gennes || can also be described by a GL model. 

For small e 2 or large Ginzburg parameter k, the GL-model lies in the deep type II 
regime where it goes over into a pure complex |^| 4 -model which describes the superfluid 
phase transition, which is the best-understood second-order transition experimentally and 
theoretically. We may therefore expect a large-/? portion of the type II regime to undergo a 
second-order transition. Indeed, the Monte Carlo simulations of Ref. |6| and the non-linear 
realization of Ref. [7] were both in this regime. 

The arguments leading to a first-order transition could therefore at best be reliable in 
the large-e 2 regime where the superconductor is of type I and the Ginzburg parameter k is 
small. Only in this regime can the e-expansion be trusted and the important problem arose 
how to find a theory explaining both type I and type II regimes. 
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One solution to this problem was proposed in 1982 by one of usfLOl, |TT| by deriving a 
disorder field theory dual to the original GL model. The RG analysis of the disorder field 



theory shows very clearly that a charged fixed point exist in the original GL model ||12|| . The 
disorder field theory exhibited a first- and a second-order regime with a tricritical point atJT0| 
n t = (3/7t)a/3/2 [1 — (4/9) 4 ] = 0.79/a/2, a result confirmed only recently with excellent 
agreement by elaborate Monte Carlo simulations |TJJ which gave K t = (0.76 ± 0M)/V2. 
Since the disorder theory is in pure |0| 4 universality class, the critical exponent v of the 
disorder field in the second-order regime has the XY- value v m 0.67 [|12[], a result confirmed 



by Monte Carlo simulations |15| . 

The RG situation of the original GL model, however, remained completely unclear. No 
fixed point at N = 2 appeared in an extension to two-loops by Tessmann a result 



confirmed later by Kolnberger and Folk using a more thorough analysis [17, IS]. Folk and 
Holovatch [|19[ managed to find a charged fixed point only by setting up a suitable Pade- 
Borel resummation of the e-expansion. A three-loop calculation using the e-expansion is 
presently under progress. As a first step, the ground state energy of the GL model was 
computed recently up to three-loop by Kastening, Kleinert, and Van den Bossche [2(J . 

Another method that leads to a charged fixed point at N = 2 was used by Bergerhof et 
al. ||21|| . Their method is based in the so called exact renormalization group ||22|| . In this 
approach, exact renormalization group equations are solved approximately using a truncated 
form for a scale-dependent effective action. The problem with this type of approach is that 
it seems to be uncontrolled, that is, there is no obvious expansion parameter. But this 
seems to be a common feature of most fixed dimension approaches. However, as we shall 
see later on in this paper, fixed dimension approaches are the best candidates to solve this 
old problem. 



At fixed dimension d = 3 there is the approach by Herbut and Tesanovic[[23|] who per- 
formed a perturbative calculation directly at the critical point. These authors obtained 
a charged fixed point at one loop and N = 2 by using two different momentum scales 
and different renormalization points for the two dimensionless couplings, / = e 2 / fi' 4 ~ d and 
g = u/fi 4 ~ d . This introduces a new parameter c = fi/fi'. They find charged fixed points 
if c > 5.16. Their approach has the problem that the parameter c must be fixed from the 



outside. They do this by demanding that the value of the ratio k = ^Jg/2f matches at 
the tricritical point with Kleinert's value[|II]] n t ~ 0.79/\/2, and found a critical exponent 
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v w 0.53. The authors used also an older Monte Carlo estimate ||24|| K t ~ 0.42/y / 2 to fix c 
which led to a v closer to the the XY value. In the light of the recent simulation work fi~4} 
obtaining the value n t ~ 0.77/\/2, the result K t — 0.42/\/2 becomes obsolete. 

There are two important universal observations made in Ref. [23|: first, if a charged 
fixed point exists, then the anomalous dimension of the vector potential is exactly given by 
7] a = 4 — d for d e (2, 4]. Second, this anomalous dimension implies the existence of [|l|, |23| 
a dimension-independent scaling relation A oc £ between the magnetic penetration depth A 
and the coherence length £. Subsequently, the AC conductivity o~{uj) was shown f25|, [26|1 to 
have the scaling behavior <j(uj) ~ £ z ~ 2 , where z is the dynamic exponent. These results 
are in contrast to the dimension- dependent scaling relations A ~ £( d_2 )/ 2 and o~{uj) ~ ^ 2 ~ d + z 
derived by D.S. Fisher et a/.|27| by neglecting thermal fluctuations of the magnetic field, 
leading to t/a = 0. 

It is important to emphasize that in Ref. ^ the two renormalization scales \x and ^ play 
a decisive role for the emergence of the charged fixed point. However, since the correlation 
functions were calculated at the critical point, there was no direct relation between fi and 
fi' and the two physical length scales of the GL model. One of this lengths, the penetration 
depth, is observable only in the ordered, since the vector potential becomes massive only 
below T c through the Anderson-Higgs mechanism. It is therefore desirable to study the RG 
flow in such a way that T c is approached from below, and this is what we shall do in this 
paper by employing a semi-perturbative approach. 

Since below T c the vector potential is a massive field, it provides us with the second 
physical mass scale of the system, and thus with the Ginzburg parameter k on which the 
RG functions depend explicitly. We shall exhibit a charged fixed point for N = 2 and 
k > l/\/2, that is, inside the type II regime. In our analysis, An = k — l/\/2 appears as 
a natural expansion parameter when T c is approached from below. This is related to the 
existence of a singular behavior at k — 1/ a/2. This singular behavior is already apparent in 
the mean field solution in a uniform external magnetic field by Abrikosov For example, 
the magnetization diverges at k = 1/ V2. To our knowledge, in the RG context this singular 
behavior in k has never been explored before. We shall show that the existence of the 
charged fixed point is related to a singular behavior at k = 1/ a/2 occuring in the vector 
potential correlation function. This singularity has the physical meaning that it represents 
the point of separation between type I and type II superconductivity. As discussed in Ref. 
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14, the value of k separating these two regimes should coincide with n t , being therefore 



lower than 1 / y/2. The approximation we consider here is not able to correct upon the value 
K = l/y/2, but higher order corrections will do. At this point it is worth to explain what 
we mean by higher order. Our method is semi-perturbative. By this we mean that we start 
by assuming that a perturbative expansion in / and g is possible. However, we parametrize 
our theory by / and k such that all RG functions depend on these parameters. This is 
possible because k 2 = m 2 /m\ = g/2f, where m = C," 1 is the Higgs mass and rriA = A -1 is 
the vector potential mass. We use m as the running scale for the RG. Thus, k arises in the 
RG functions from two different sources: the coupling g, which is rewritten as g = 2Jk 2 , 
and the loop integrals, from where the masses are combined in such a way as to produce 
functions of k. The RG functions obtained in this way are polynomials in / with coefficients 
depending on k. Thus, it is not really a standard perturbative series. This procedure will be 
explained in Sections III and IV. The key point is that for k sufficiently close to k = 1/ y/2 
from above, / is small enough such that a charged fixed point can be obtained at N = 2. 

An interesting point revealed by our analysis is the role of k in controlling the appearance 
of momentum space instabilities in the correlation functions. For k > 1/ y/2, the two-point 
bare correlation function is maximized for a nonzero wave vector, similar to the roton peak 
in the correlation function of superfluid helium. This behavior seems to be the origin of the 
negative sign of the ^-exponent and implies a critical behavior |2B[] known from the theory of 
Lifshitz points. In contrast to scalar models of Lifshitz point, where the momentum space 
instability is already present at the tree level, the instability in the GL model is induced by 
magnetic fluctuations. Interestingly, such fluctuation-induced Lifshitz points also occur in 



the non-commutative -theory, which also possesses a negative r\ exponent [28 



The plan of the paper is the following. In Section II we discuss the advantages of a 
RG approach to the GL model in d = 3 dimensions in the disordered phase over the usual 
4 — e-dimensional calculations. This approach is a generalization to the GL model of Parisi's 
three-dimensional RG technique P9| . In particular, we show that this approach, being 
sensitive to the infrared divergences of the vector potential, yields at the one-loop order 
a more reliable result than the one-loop e-expansion. In Sections III we renormalize the 
ordered phase followed by a detailed renormalization group analysis in Section IV which 
culminates in the desired infrared- stable fixed point which could not be found in 4 — e 
dimensions. Some useful integrals used in the paper compiled in the Appendix. 
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II. THREE-DIMENSIONAL RENORMALIZATION GROUP IN GL MODEL 
AND ITS ADVANTAGES 



The bare GL hamiltonian considered in this paper is given by 

H = i(V x A ) 2 + |(V - ie A )4j \ 2 + m 2 \^ \ 2 + ^ | 4 + H#, (1) 

where 

^gf = ^(V-A ) 2 (2) 
fixes the gauge. It will be convenient to express the complex fields in terms of real fields as 

^ = ^(tf ) + z4 2) )- (3) 
The propagators are given in momentum space by 

(p) = (4\p)4 j) (-P)) = -J^-s , (4) 



p z + m, 



o 



z>jg( P ) = «(p)^(- P )) = i + (a - l)^ 



p 2 



(5) 



Since the free vector field is massless, perturbative calculation run into infrared problems. 
For instance, the Feynman graph in Fig. 1 contributing to the four-point function is infrared 
divergent at zero external momenta for any dimension d G (2, 4). Indeed, the loop integral in 
this graph is proportional to that the Feynman integral yields in dimensional 

regularization |I7|, |3D| a pole term 1/e (e = 4 — d). However, this pole reflects an 
ultraviolet divergence for d — 4, not an infrared divergence, which is physical source of 
critical properties. Recall that the reason why ultraviolet divergences give nevertheless 
information on the infrared behavior is that a 1/e pole in d = 4 — e is equivalent to a 
logarithmic divergence of the form ln(p 2 /A 2 ) at d = 4, where A is the ultraviolet cutoff, and 
the ultraviolet limit A —>■ oo at fixed |p| is equivalent to taking the infrared limit |p| — >• at 
fixed A. 

Alternatively, we can choose |p| = /i ^ as a scale parameter of the problem, confining 
the analysis to d G (2, 4) to avoid the logarithmic divergence at d — 4. This procedure is 
applicable in a RG analysis at the critical point performed by Herbut and Tesanovic [p3| . 
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However, as we have discussed in the introduction, this procedure needs a second scale 
parameter // which has to be determined from the outside. 

Another way of regulating the infrared divergences of the GL model is to introduce a 
Proca term in the Hamiltonian, Mq Aq/2, which explicitly breaks gauge invariance. Then 
we may perform the calculations of correlation functions using the infrared-finite propagator 



£J2(p) 



p 2 + Ml 



8uv + (oc - 1); 



(6) 



"" v u V + a Q M 2 
and taking Mq — > at the end. However, if dimensional regularization is used to evaluate the 

massive Feynman integrals, final results are the same as in the dimensional regularization 

of the massless case, since the 1/e pole is the same for all M : 

d d p 1 1 



(2vr) d (p 2 + Ml 



2\2 C^2 



+ 0(1). (7) 



Thus, the /3-functions obtained with this procedure will have the same e-expansion as in the 
original work||, where a charged fixed point exists only for a number of ip^ field components 
N > 365. 

A completely different result is obtained if we study the system directly in three dimen- 
sions. We specify the renormalization constants by rewriting the Hamiltonian (fl|) in terms 
of renormalized quantities as 

H = ^(V x A) 2 + ^A 2 + Z„|(V - *eA)Vf + Z m m a \i/>\ 2 + Z g ^\ A + H gi , (8) 

where the renormalized fields are defined by tp = tjj and A = Z A A . We shall work 
in the Landau gauge where ao = in 7i g f. The term eA does not need a renormalization 
as a consequence of the Ward identity |T| which implies that e 2 = Za^- The renormaliza- 
tion constants Z^, Za, Z m , and Z g and the renormalized Proca mass M are fixed using the 
normalization conditions for the one-particle irreducible renormalized n-point functions [|I|, [2j 
(no sum over repeated i indices): 

dii 2) ( P ) , 



lp=0 



1, (9) 



dp 2 

ri 2) (o) = m 2 , (10) 

ri(0, 0,0,0) = 3ms, (11) 
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ar$(p) 



<9p 2 



(12) 



p=0 



rS(o) = 3M 2 . 



It is straightforward to calculate the following results at the one-loop level: 



1 + 



3vr m + M ' 



(13) 



(14) 



Z A 



24vrm' 



(15) 



5w 2 



m 3irm(m + M) Sum 2 2nmM 
We define a dimensionless renormalized square charge by / = e 2 /m. Using e 2 
find the one-loop equation 

p 2 4 

j — _2 _ 
m 247rm 2 

The /3-functions are defined in general by the logarithmic derivatives 

* - r df 
at = hm m— — , 

m->o dm 



(16) 
Za^q we 



(17) 



(18) 



/3„ = lim m 



dg 



(19) 



Af->o 9m 

The derivatives with respect to m have to be performed at fixed bare couplings. At this 
point we make a crucial observation. In the derivatives, we reexpress the bare couplings in 
terms of the renormalization / and g, while discarding terms beyond one-loop. However, the 
term proportional to in Eq. (ITH) contains only a single power of m in the denominator. 
Therefore this term will contribute only to the term —g in (3 g and no term proportional to 
f 2 will be present. The result is 



rn 



dg 
dm 



u 



elu 



+ 



5u 2 



+ 



m 3ir m(m + M) 3tt (m + M) 2 Airm 2 2irmM 



elu 



m 3nm(m + M) 
4 



5m 2 



mm' 



2-nmM 



elu 



+ 



5m 2 



3tt (m + M) 2 Mm 2 



-9 



e 2 u Q 



+ 



5u 2 



3tt (m + Mf 8nm 2 



(20) 



The derivatives of M with respect to m do not contribute since they are of higher than 
one- loop order. Taking the physically relevant limit M — > 0, we obtain 

Pg ~ y o 9 ~* 



3n m 2 %7im 2 

4 . V 

671 8/T 

where in the second line we have replaced u /m —>■ g and e\jm — > f, the error committed 



4 5o 2 

-g-Trfg + T, (21) 



with this substitution is beyond the one-loop order under consideration. It is this result 
where the three-dimensional calculation in this paper differs essentially from the dimension- 
ally regularized one in d = 4 — e dimensions in which the /3-function j3 g contains additional 
f 2 term. Precisely this term is the culprit for the nonexistence of a charged fixed point in 
HLM. 

The second /3-function is 

/>/ = -/ + £. (») 

From Eqs. (|2l|) and fl22[) we obtain the infrared stable fixed point /* = 2Ati and g>* = 2647r/5. 

At this point we may think that the problem is solved, and that the absence of a charged 
fixed point was merely an artifact of the e-expansion. Unfortunately this is not true since at 
the fixed point, the critical exponents have completely unphysical values. For example, the 
critical exponent rj, which is defined by the fixed point value of the RG function 

7,7, = lim ra ^ln^ 

v m^o dm 

has the value 77 = — 16, which is unphysical since it does not respect the bound 77 > —1. The 
reason for this is that the fixed-point value /* is quite large, and this leads to the unphysical 
value of 77. 

Does this mean that the present approach must be discarded as well? Not completely. 
We are still much better off than HLM. As mentioned before, they judged how far their 
result was from nature, by changing the number of complex fields from 1 to N/2 with an 
O(N) symmetry among the real components. This served to suppress the f 2 term, since the 
graph of Fig. 1 is of order|36| 1/N 2 . Our calculation has the advantage that the f 2 term is 
absent for any N. Suppose we go through our analysis for N > 2. Then we find 

Nf 2 

ft = -' + lSF' (24) 
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Afg (N + S)g 2 

^ = -^-^T + ^6^- (25) 

The critical exponent has now the value i] = — 32/iV, which is physically meaningful provided 
N > N c = 32. Thus our three-dimensional approach yields a critical N value N c = 32, where 
the result becomes physical. This lies much lower than the critical HLM one-loop value |3j 
jyHLM = 365. If our one-loop value N c = 32 is lowered by the two-loop corrections, there is 
a good chance of getting a physical exponent rj for a single complex order field ip. Such a 
calculation remains to be done. 

Further insight can be obtained by calculating the critical exponent v. This is given by 

1 = 2 + 7^-7?, (26) 
where is the fixed point value of the RG function 

= m 77~~ • (27) 



At one-loop order it is given by 



dm 



Thus 



1 (^ + 2)^ + 64) 16 

v N(N + 8) N K ' 

Positive values of v are obtained for iV > 34. However, for any finite value of N we have 
v > 1. When iV — > oo we have v — 1. 

As a last remark in this section, let us comment about the absence of a tricritical fixed 
point in our one-loop calculation. The tricritical fixed point is absent because there is no 
f 2 in the /3-functions ( plf ) and (p5|) . With this respect, our approach gives a result similar 
to the 1/N expansion at 0(1/ N), where there is only a second-order fixed point and no 
tricritical fixed point. It is evident, however, that a tricritical fixed point will be generated 
at two-loops, but in this case resummation is necessary. The same result is to be expected 
in the 1/N expansion at 0(1/N 2 ). 
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III. RENORMALIZATION CONSTANTS IN THE ORDERED PHASE 



A. Problems with unitary gauge 

We have seen in the preceding section that in three dimensions the situation improves by 
using a massive vector field to avoid infrared divergences in the Feynman integrals. Above T c , 
this mass is set equal to zero at the end. In the ordered phase, such a mass is automatically 
present as a consequence of the Meissner effect. In fact, there are no massless modes at 
all in the ordered phase, the Goldstone boson supplying the longitudinal component to the 
massive vector field. The simplest way of representing this effect is in the unitary gauge. At 
the mean-field level, this amounts to the field parametrizations (dropping the subscripts 
indicating bare quantities for simplicity) 

^ = -L(v + p)e w J (30) 

B = A + -V9. (31) 

e 

The above parametrization allows for a nonzero expectation value of the order field, that 
is, by setting (-0) = v. At the tree level this happens for m 2 < 0. To avoid proliferating 
minus signs, it will be more convenient in this section to replace the term m 2 \%l)\ 2 in the 
Hamiltonian by —m 2 \tp\ 2 , such that a nonzero expectation value of the order field exists for 
m 2 > 0. In the Coulomb gauge V ■ A = the partition function is calculated from the 
following functional integral 

Z = J VAVpVOp det(-d 2 )5(V ■ A)e-f d3rn , (32) 

the last factor in the measure being the Faddeev- Popov determinant. This is a convenient 
formulation of the field system since by performing the change of variables ( |30"D and (|3T| ) the 
functional integral can immediately be simplified by integrating out the angular field 9. Its 
value is fixed by the delta function enforcing the Coulomb gauge. The result of this integral 
precisely cancel the Faddeev- Popov determinant. 

In terms of the fields (|30"D , the Hamiltonian becomes 

U = V (v) + ~(V x B) 2 + ^B 2 + \iVpf + ^p 2 B 2 + e 2 vpB 2 + m 2 p 2 + ^p 3 + ^p\ (33) 
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where Vq(v) = —m 2 v 2 /2 + uv 4 /8, m\ = e 2 v 2 , and m 2 = uv 2 /2. In this well known tree level 
analysis the Goldstone has disappeared and the vector potential has becomes massive. 
The free propagator of the vector field is given by 

^(P) = z5-^rf^ + ^y ( 34 ) 



which corresponds to take the gauge-fixing parameter a — > oo. Now, at large |p| the above 
propagator behaves like a constant. Thus, the free scaling dimension of B does not have 
the typical canonical scaling dimension 1/2 of a vector field, but 3/2. But then power 
counting tells us that the interaction e 2 p 2 B 2 /2 is not renormalizable in d = 3. Thus, in the 
unitary gauge, the absorption of the Goldstone mode in the longitudinal part of the massive 
vector meson seems to destroy renormalizability. The only way to obtain a renormalizable 
theory out of the unitary gauge is by keeping the gauge-fixing parameter a nonzero during 



the Feynman integrals calculation and take the limit a — > oo at the end [54]. In order 
to implement such a program, we should not use the unitary gauge parametrization given 
by Eqs. (|30|) and flST]). This means that the Goldstone boson will still be present in the 
theory. The pole coming from the Goldstone boson propagator will be, however, canceled 
from physical quantities. 

Beside the renormalizability problem, there is another difficulty with the unitary gauge 
which is related to the vortex content of the theory. The point is that the covariant derivative 
term |(V — ieA)ip\ 2 in the GL Hamiltonian (|I|) goes over, in the parametrizations fl3!]| ) and 
(j31~l), into |(V — ieQ — 2ttO v )p\ 2 , where V is the vortex gauge field, a vector field describing 



the vortices in the superfluid [|IT], [13]]. This field arises from the fact that 9 is not a single- 
valued field, but a a cyclic field with the property J dO(r) = 27rn„(r), where n v is the winding 
number or vorticity. For this reason, the gradient Ve ie is not simply equal to i(V9)e ld , but 
must be supplemented by such a vortex gauge field yielding i{V9 — 2n6 v )e' td , and ensuring 
the periodicity under 9 — > 9 + 2tt. The vector gauge field removes possible delta functions 
arising from the gradient of jumps by 27r in the cyclic field which exist around a vortex line 
in the superconducting phase [JXTJ, [13]]. In fact, these vortices are crucial for driving the phase 
transition and can therefore not be neglected. It is, however, difficult to treat the vortex 
gauge field by conventional perturbation theory using Feynman diagrammatic techniques. 
If we want to take vortices into account in perturbation expansions we must therefore avoid 
multivalued fields. 
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B. The Landau gauge 



A field parametrization which avoids the problems of the unitary gauge and has more 
desirable properties with respect to power counting is 

tpo = -j=(v + o- + in ). (35) 
V2 

Here the bare Hamiltonian becomes 

_ e 2 

ri = ^frcc + e A • (o-oVtto - 7T Vo- ) + e 2 v a Al + -^(cr 2 + 7r 2 )A 2 

U V 3 U V 2 , U 0f 2 , 2^2 foa\ 

+ ^-^o + —°on + -g-(<7 + tt ) > ( 36 ) 
where 'Hfree denotes the free part of the Hamiltonian: 



' V X A ) 2 1 "' Kil A 2 1 ' "TT_ \2 , /T7„ \2l , ' , o.,.2 ,„2 

1 



H frcc = -(V x A ) 2 + -^A 2 + -[(Vao) 2 + (Vvro) 2 ] + -(-m 2 + 3m 2 )a 2 



+ ^( _m o + '"o)^ + J o°o + ^ g f- (37) 

The Coulomb gauge V- A = is fixed by letting a - ► 0, which corresponds to the so called 
Landau gauge. We have therefore omitted a crossed term e A • V7r . The bare masses are 
m o — UqVq/2 and m 2 A0 = e^v 2 . We have introduced a source term for the longitudinal field 
ctq. In the equation of motion, a zero source corresponds to the minimum of the effective 
action. At the tree level we have Jo = vo(—ml +Wq). Thus, the tree level minimum implies 

TTIq = TTIq. 

From Tifree we obtain the propagators: 

3S(P) = 2 * 2 _ 2 , (38) 

p 2 + 3mg — TTIq 

G<g(p) = 2j _ \ _ 2 , (39) 
p 2 + - 



^«P) = -tA^" f ^ " ^ ) ■ (40) 



^ p 2 + mi i0 V^ P 2 
Now power counting gives the desired dimensions in the ultraviolet since the vector prop- 
agator behaves like ~ 1/p 2 for large |p|. Moreover, the vector field is massive and the 
graph of Fig. 1 is convergent in d = 3. However, we still have a massless mode, the 7r-field. 

13 



Interestingly, in the one-loop calculation of the vector potential and ix two-point functions, 

(2) (2) 

Tfj,v and T^, respectively, the infrared divergences from the would-be Goldstone boson do 
not play any role. We shall calculate here the bare two-point functions. The Feynman 
graphs contributing to these two-point functions are shown in Fig. 2. The only ultraviolet 
divergences come from the tadpoles and they are all proportional to the ultraviolet cutoff 
A. These divergences are absorbed into the bare masses Wq and m\ . Also, in the loop 
integrals we shall replace the bare masses by the renormalized ones, since the error involved 
in this replacement is of higher order. 

The renormalized two-point functions are given in terms of the bare ones by T^l = 
^a^ojaai = ^r^o^Tr) anc ^ ^jS = Z 'aTqui/- in order to determine the renormalization 
constants Z v and we employ the normalization conditions 

= 1, (41) 



<9r$( P ) 



dp 2 

ar$(p) 



p=0 



= 2. (42) 

p=0 



dp 2 

Since 7r is a would-be Goldstone boson, we require also the normalization conditions for 
renormalized source J = Z^ 2 J = 0: 



r^(0) = uv 2 = 2m 2 , (43) 

rg}(0) = 0, (44) 

such that all 7i~-propagators in loop integrals are massless. 

The following normalization condition defines the renormalized mass of the vector field 

rg(0) = 3cV = 3< (45) 

Note that, as usual, the renormalization conditions are chosen in such a way as to be con- 
sistent with the tree-level and the Ward identities [0. Note that the above conditions 
determine completely the renormalization constants. Since we have five renormalization 
constants, the above five conditions determine them completely. Actually, the Ward identi- 
ties imply that Z a = Z n if the Coulomb gauge is chosen, a result consistent with the gauge 
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invariance of the Hamiltonian. The renormalization conditions imply that the renormalized 
Ginzburg parameter is given by k 2 = m 2 /m 2 A = u/2e 2 , that is, it has the same form as at 
the tree-level but with the bare couplings replaced by the renormalized ones. 
After some work we obtain 

r (2) , s _ p 2 ^2 . m 2 u V2m e 2 m A 



u m 2 
47r|p| 



6(p, V2m) - el 1 - — \—-[Am 2 A p 2 + (p 2 + 2m 2 - m^) 2 ]a(p, V2m, m A ) 
y&ixm A \p\ 



+ ilb-ra (p2+2m2)2Kp ^ m) ^ (46) 



Op) = V \ p 2 + <„ - ^(V2 K - 1) - «(i - 2^ - V2 K) e »^'" 3 



87r 87rp 2 47rp 2 



+ 



e 2 [p 2 + rrfl( V2« + l) 2 ] [p 2 + m\ (V2k - l) 2 ] - Smjp 2 /- ' 
16^ p a(p ' V2m ' mA) 



+ i^{_ p2 + ^ K _3 ) + M^! + |M (1 _ 2K2 _^) 



2 4 2 

e^m , , s ef 



6(p, v^m) - — ^-[3(p 2 + 2m 2 ) 2 - 12m 2 m 2 4 



27r|p| 3 ' 16/T | p| 3 

+ 3m A -2m 2 A p 2 ]a(p,V2m,m A )} (47) 

with the functions 

/ x fp 2 + m 2 — ml\ fp 2 + ml — m 2 \ 

a(p, mi,TO 2 J = arctan : — : + arctan : — : , (48) 

V 2m 2 |p| ) \ 2mi|p| / 

b(p, mi) = lim a(p, m 1 ,m 2 ) 

/ 2 2 
7T IP m l 



+ arctan , (49) 

2 V 2mi p - 
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where k = m/m^ = A/£ is the renormalized Ginzburg parameter. The bare Ginzburg 
parameter is given by k = a/ u /2el, since ml = UqVq/2 and m\ A = e^v^. 

On account of the normalization condition fl44|), we obtain the one-loop correction to m\. 



m 2 =m 2 - 3UoV * m - (50) 

87T 27T 

In the following, we shall use extensively the replacement Kq — > k 2 = m 2 jm\ = u 2 /2e 2 = 
g/2f, where f = e 2 /m and g = u/m are dimensionless renormalized couplings, neglecting all 
errors of higher than one-loop order. The parameter k arises in two different contexts: once 
from the renormalized mass ratio m/m^, and once from the ratio of couplings as \Ju/2e 2 . 

We shall then parametrize all RG functions in terms of k and / instead of g and /. 
Our expansion will be controlled initially by powers of /, but once g is eliminated in favor 
of k, then k will no longer be assumed to be small, such that the RG functions are not 
polynomials in / and k as in the previous expansions in powers of / and g in Section II. 

Below we shall see that in this formulation there exists a natural expansion parameter 
related to k, namely Ak = ac — 1/ v2- Remarkably, the present description of the fixed-point 
structure requires only the knowledge of the two-point functions, since from Eq. ([50]) m 2 = 
Z^rh 2 and e 2 = Za&q, with Za being determined from the renormalization condition (fffi) (see 
below). This represents an immense advantage since the four-point functions contain severe 
infrared singularities coming from the would-be Goldstone boson. Of course, we expect these 
singularities to cancel at the end, but this is a difficult issue in fixed dimensions. In the 4 — e 
-dimensional regularization scheme this issue is avoided since the counterterms required for 
a renormalization in the ordered phase are identical to those in the disordered phase where 
they can be calculated without Goldstone-bosons [![]. 

From the normalization condition ( f4"T| ) we derive 

12tt( v / 2k + 1) 

which satisfies Z^ > 1 for k < k = (V33 - l)/2y/2 w 1.6774. Only for k > k will the wave 
function renormalization satisfy the usual inequality < Z^ < 1 found in the absence of 
gauge fields. 

Let us now calculate m\ at one-loop level using the normalization condition (^5J). We 
find 
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Z?mW (52) 



Since T$ = Z^T^ , we obtain 



m A = Z A Z n m AQ . (53) 

Thus Z n appears in the renormalization of the mass m A of the vector field. This result is 
expected from the Ward identities. It provides us with a good check of the consistency of 
our calculations. Physically it shows clearly how the fluctuations of the Goldstone boson 
renormalize the mass of the vector field, which has been exploited in the theory of electroweak 
interactions to build a renormalizable theory of massive vector mesons and Higgs fields. 
Thus, although the Goldstone boson is present in the calculations, the Higgs mechanism 
takes care of absorbing its fluctuating degrees of freedom to build a fluctuating longitudinal 
degree of freedom for the vector potential. 

The superfluid density is intimately related to the penetration depth A = l/m^. It is 
given by 

Ps = Z-\l (54) 
which reflects a relation due to Josephson[]3~8fl . On account of Q53|) and the relation e 2 = Z^e^ 



we can write m\ = e 2 p s . 

It remains to calculate Za- This is done using the normalization condition (fE|), yielding 

- - i ^ C(K)/ (55) 



24vr(2/t 2 - l) 3 
where 



C(k) = 4k 6 + 10k 4 - 24V2V + 27k 2 + AV2k - 1/2. (56) 

We note that the second term in Eq. (|55]) is singular at k — 1/ y/2. This singularity is 
analogous to the 1/e singularity in dimensionally regularized theories. The important role 
of this singularity will be discussed in more detail in the next Section. 
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IV. RENORMALIZATION GROUP ANALYSIS IN THE ORDERED PHASE. 
ROLE OF THE GINZBURG PARAMETER k 



A. Renormalization group functions 

We are now prepared to calculate the RG functions. Let us define the RG functions: 

7 ^ = m ^T' (57) 

_ dlnZ A 

7A = m^—. (58) 
The RG functions 7^ and 7,4 are given explicitly at one-loop order as 

^ = iSF ^ + 1 ' (59) 

7 ^ = 24^(2^-1)3 - ( 60 ) 



Note that while deriving Eqs. (|59| ) and ((6^) we have not differentiated k, since this would 
lead to higher powers of / beyond the one-loop approximation. We observe also that the 
singularity at k = l/y/2 present at Z A is also present in 7^4. 
The /9-function of / is then given by 

P f = (ta - I)/- (61) 

To obtain the /3-function of k 2 we need to know the evolution equation for m\. From Eq. ( |53|) 
we obtain 

thus reducing the problem to a calculation of mdm\ /dm. It must be emphasized that in 
our approach this is not zero. Our differentiations are performed at fixed Kq = m 2 ) /m 2 i0 , 
such that mdnl/dm = whereas both m and m^o remain functions of m. Due to the 
fixed Kq, their derivatives mdm 2 A0 /dm = (m 2 i0 and mdm^/dm = C m o> are governed by a 
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common RG function (. In order to obtain mdm^j 'dm we use Eq. ([50]) to write 



m 2 _ z -i 2 , Wjjm e 2 m A 
m -Z n m + ^ + -^ r , (63) 



where we have inserted m 2 = Z n m 2 . The derivative yields 

dm 2 



m- 



dm 



^ = (2 + Cr-TirK, (64) 



where up to one loop 



4vr V 2 v 7 ^ 

From Eqs. ( |52] ) and (p5 ) we derive 



C--^/(^ + ^). (65) 



= (2 + C + 7a - 2 77r )m A . (66) 
It is now straightforward to obtain the /3-function for k 2 : 

8k 2 dm 2 /m 2 A ( m dm 2 A \ , . , . 

f3 K 2 = m— = m J A =k 2 [2- = k 2 (2 1it -<y A - Q. (67) 

am am \ m A am J 

From Eqs. (|6l|) and (|67D we see that charged fixed points satisfy the equations 

Va = 7a(/*, «*) = 1 (68) 

and 

277 - 1 = Ctt (/*,«*), where 77 = 7^(7*, /«*)• (69) 

The exponents 77^ and 77 determine the anomalous scaling dimensions of the fields A and n, 
respectively. Our equations yield fixed points at 

/♦»0.3, k,«1.17/V2. (70) 

At this fixed point the value of the 77-exponent is 77 ~ —0.02, which fulfills the inequality 
7] > — 1, showing that our fixed point in the ordered phase is completely physical, in contrast 
to the N = 2 case in the disordered phase calculation of Section II. 

Our ?7-exponent is less negative than most values found in the literature. There is, 
however, no consensus about this value. For instance, a recent Monte Carlo simulation 
3"4| give 77 ~ —0.24, disagreeing with another Monte Carlo value 77 « —0.79 published 



19 



recently |35[], which is similar to i] = —0.74 found by Herbut and Tesanovic by adjusting their 
parameter c to fit the Kleinert's value of the tricritical Ginzburg parameter, n t = 0.77/v^2- 
The value rj = —0.2 is obtained when the c-parameter is fitted to the now less precise (see 
Ref. 0) Monte Carlo value [|2§ n t = 0.42/ v^. There are other values of rj reported in the 



literature JTL], 33[], exhibiting a lack of concensus on its numerical value. 



B. Validity of the method 

We now want to explain what makes our three-dimensional procedure so successful in 
finding the charged fixed points. The main obstacle in our calculation in the disordered 
phase in Section II was the too large value of /*. In previous approaches in 4 — e dimensions, 
this prevented a fixed point to exist. The size of /* was diminished by a factor 2/N via 
the artifical extension to N/2 complex fields, and this led to a fixed point after all for 
iV > N c = 365. The three-dimensional calculation in the disordered phase is far less 
unphysical since it always give a fixed point. Still, there was the problem that the associated 
critical exponent r\ was too negative, violating the bound r\ > — 1, again because of a too 
large /*, unless N > N c = 32. 

These difficulties are absent in three dimensions below T c . To understand this let us 
define an effective square charge f(n) depending on k by the equation 

7a(/,«) = 1, (71) 

meaning that we go to the "fixed line" in the two-dimensional space of coupling constants 
defined by the vanishing of the /3-function /3f in Eq. (|H]), but not at the fixed point for Eq. 
( |S7|) . The ^-dependence of /(/«) is 

= 24.(2.3 -I)* 
V ; V2C{k) V ; 

We see that f(n) = for k = 1/V2, which is just the value separating the type I from 

the type II regime at the mean-field level. The ^-dependence of J(k) is plotted in Fig. || 

The effective square charge is neg ative if 0.096/^ < k < l/y/2 and very large positive for 

< k < 0.096/v2- This makes it impossible to have an infrared-stable positive charge in 

the type-I regime k < 1/ y/2, and we conclude that the phase transition is of first-order in 
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the type-I regime. Moreover, we observe that for k slightly above l/\/2, the effective critical 
charge is small. 

An important remark is in order here. In Fig. |3| there is an asymptote which separates 
two distinct regions: one where the charged fixed point is inaccessible (the left side of 
the asymptote), and another one where the charge fixed point is accessible. This asymptote 
reveals a Landau-ghost-like behavior. The so called "runaway" of the RG flow[|3| corresponds 
to the left side of the asymptote. It is worth mentioning that a similar feature can be 



expected in the /3-function of QCDfl39|. There the asymptote is supposed to separate the 
asymptotically free regime from the regime containing an infrared stable fixed point which 
governs quark confinement. 

Consider now the large-K limit which corresponds to an extreme type II regime. There 
are two ways of taking k — > oo: either by letting g — > oo at / fixed, or by letting / — ► 
at g fixed. In the first case we obtain from Eq. ( (58]) , Ja\k=ooJ — //(247rv / 2), while in the 
second: 7aU=oo,p = 0. In the first case the magnetic fluctuations are still relevant and the 
\ip | 4 interaction may be replaced by a pure phase field with fixed absolute value of the order 
parameter, i.e., the GL model be be approximated by its London limit. In the second case 
the magnetic fluctuations disappear completely from the system and the universality class 
is that of a pure |?/>| 4 -theory which is equal to that of the XY -model. 

By substituting /(/«) of Eq. ([72] ) into (^9|), we obtain the effective exponent For 
k = l/V%, this has the mean-field value 77 = fj(l/^/2) = 0. In the physically interesting 
interval l/v^ < k < k = (V33 — 1)/2a/2 w 1.67746, the effective exponent fj(K) is negative, 
as shown in Fig. f|. Above k, it is positive. Remarkably, it lies above the physical lower 
bound —1 for all k > 1/a/2, i.e., in the entire type-II regime in the mean-field approximation. 

Thus, our RG approach in the ordered phase gives effectively an expansion around k = 
1/V2. It is therefore convenient to introduce the expansion parameter Ak = k — l/v^2- 
Then we can write the leading expansion term in Ak for the effective square charge f(n) 
and the effective exponent t/(k) as 

f(n) = 487tAk 3 , (73) 
fj{K) = -6V2Ak 3 . (74) 
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Similarly, we expand 



(3 k 2(k 2 , /(«)) = -- - V2Ak - An 2 + 162v / 2Ak 3 . (75) 

We now discuss the other important critical exponent v. Recall that it is defined in 
general by the scaling relation 



v 



2 f]m 

where rj m is the fixed-point value of the logarithmic derivative 



(76) 



wig am 



the right-hand side following from fl64j) . By analogy with the other effective quantities we 
define the effective exponent 

PW = 2 + C.W -,(«)' (78) 

where Ct(k) = Ctt(/( k ))- The effective exponent P(ac) gives the critical exponent z/ at the 
fixed point k = ft*. To leading order in Ak we obtain 

P(«) = \ + ^|A. 3 . (79) 

By substituting the previously found fixed-point value = 1.17/v^2 of Eq. (|70|) into 
([78|), we obtain v pa 1.02. To estimate the systematic error, we calculate the fixed-point 
value k* from the expanded /3-function ([75]), which yields = 1.22/V2. Inserting this into 
Eq. flTJ) gives i/ w 0.92. The expanded Eq. (|7|) yields 77 pa -0.03. 

The critical values of v are quite far from the expected XY-model value vxy ~ 0.67, but 
at the one-loop level we should not expect a higher accuracy. 

Systematic improvements in our approach is considerably more difficult than in other 
more conventional approaches. The fundamental difference between our method and other 
methods lies in the fact that we perturbatively expand in powers of / only. The essence of 
the method lies in the computation of and further determination of / using Eq. (f71~l). 
Generally we have that 7^ is given by the series: 

lA (f ) K) = Y,^)f. (80) 
1=1 
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The powers of / correspond to the number of loops and the coefficients q(k) are not polyno- 
mials in k and should diverge at some value k = n t separating type I from type II regimes. 
At one-loop order n t = 1/ v2, receiving no corrections with respect to the mean-field value. 



C. Scaling behavior as T c is approached from below 

Let us discuss the scaling behavior as T c is approached from below from a perspective 
independent of perturbation theory. It is convenient to write the formulas for any dimension 
d G (2,4). Thus, the /3-funtion for the gauge coupling is given by 



(3 f =(jA + d-4)f. (81) 
The /^-function j3 K 2 can be written as 

&» = ^-d- lA + ^jn 2 . (82) 

Since m\ = m 2 /n 2 , we obtain 

The existence of a charged fixed point implies tja = l\ = 4 — d an d = 0- Thus, in the 
neighborhood of this charged fixed point Eq. (|83|) becomes mdm 2 A / dm ~ 2m 2 4 , implying 
that m 2 A ~ m 2 . Therefore, A and £ diverge with the same critical exponent ]I5L [23j. This 
is of course obvious from the definition of the Ginzburg parameter: k 2 = g/2f = A 2 /£ 2 . If 
g* and /* are different from zero the ratio A/£ is a constant at the critical point. This can 
only happens if they have the same exponent. At the XY fixed point, however, tja = 
because /* = and mdm 2 A / dm ~ {d — 2)m' A near the fixed point. In this case m 2 A ~ m d ~ 2 
and the penetration depth exponent is given by v' = v[d — 2)/2. This is the so called 



XY superconducting universality class |27]. The universality class governed by the charged 
fixed point is fundamentally different from the XY universality class. Note that the above 
argument works only if we approach T c from below, since m 2 A ^ only for T < T c . 

From Eqs. ([53"D and (|54|) we can write m 2 A = e 2 p s . Due to Eq. (|8TD, e 2 ~ m 4 ~ d near 
the critical point governed by the charged fixed point. Since there m 2 A ~ m 2 , we obtain the 



Josephson relation [p8] p s ~ m 
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V. MOMENTUM SPACE INSTABILITIES. ORIGIN OF THE NEGATIVE SIGN 
OF THE I/-EXPONENT 



A much debated property of the superconducting phase transition is the sign of the 
^-exponent fZB, I3T], BBI, p 



J. Recently, it was pointed out by one of us[|26| (F.S.N.) 
that the origin of the negative sign of 77 lies on momentum space instabilities of the order- 
parameter correlation function arising from magnetic field fluctuations. Such momentum 



space instabilities are very similar to those occuring in scalar models of Lifshitz points [40 



In these models, the ^-exponent (which in this context is usually denoted by 7^4) is also 
found to be negative. The only difference with respect to the GL model is that in the scalar 
models for Lifshitz points the momentum space instabilities are included from the beginning 
explicitly into the Hamiltonian. 

,(2) 



Inspired by works on lattice spin models of Lifshitz pointsjll]], we expand up to 

order (p 2 ) 3 : 



,(2) 

0,7T7T 



r-1 



+ 



Af p 2 , Bf ( p 



+ 



7T m 



7T 



m, 



(84) 



where 



.4 



1 - V2k)(V2k + I)- 3 ( \ + IV2k + 'pi , 

\o 5 5 / 



B 



1 



l + 5 v / 2K + 20fi; 2 + 18v / 2«: 3 



(85) 



(86) 



56^ 105^ (v^^+l) 5 
The above expansion is like a Landau expansion where p play the role of an "order param- 
eter". Remarkably, the coefficient A vanishes for k = l/y/2. In the type I regime A has a 
positive sign while in the type II regime A is negative. The coefficient B is always positive, 
ensuring the stability. Close to k — l/y/2 we find 



1 ( 1 

A ~ k = 

24 V y/2 



An 
HI' 



(87) 



B 



96^' 



(88) 
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For k £ l/y/2, T^ n7T (p) has a zero at 

|po| ~ 2 3/4 m A A/^, (89) 

with the exponent (3 = 1/2. Note that r^^p) is just the inverse of the transverse suscepti- 
bility Xt(p)- Thus, the transverse susceptibility is maximized at | p» 1 - This is the behavior 



assumed so far in all generic models exhibiting a Lifshitz point |4(J. It implies the existence 
of a modulated regime for the order parameter. On the basis of these considerations we 
conclude that the transition from type I to type II behavior happens at a Lifshitz point. 

This conclusion is, in fact, in accordance with another physical property of the transition 
from type I to type II superconductors. In the second, fluctuating vortex lines play an 
important role. This is precisely the reason why the fluctuation term ~ |^| 4_e calculated 
under the assumption of a constant |?/>| is unreliable. At the core of each vortex line, \if)\ has 
to vanish, such that in the type II regime the order field os perforated by lines of zeros |x| . 

In our one-loop calculation the k- value where a type I superconductor crosses over to type 
II has the mean-field values Kl = l/y/2. This value is expected to decrease in higher- loop 
calculations towards the tricritical value K t ~ 0.77/v^2 where we expect the cross over to 
take place. Indeed, in the dual theory [K| pj| , the tricritical point comes about by the 



sign change of the quartic term in the disorder field which corresponds to a change of the 
average interaction between vortex lines changing from repulsion to attraction. This sign is 
also what distinguishes the type I from the type II regime experimentally. 

It remains to be checked by going to higher loops, that the fluctuation corrected value 
of does indeed coincide with the Ginzburg parameter at the tricritical point, that is, 
k l -0.77/^ MM- 



The exact value of the exponent (3 controlling the vanishing of |po| at the type- 1- type II 
boundary in Eq. ([59]) is expected to increase slightly above 1/2. 

An interesting open problem is to prove the conjecture of Ref. that the Lifshitz point 
coincides with the tricritical point obtained from the disorder field theory HI]]. 



VI. CONCLUSION 



In this paper we have initiated a new RG approach to the GL model. Since we work 
directly in d = 3 dimensions and in the ordered phase, we were able to express the RG 
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functions as functions of the Ginzburg parameter k = m/rriA- This replaces the coupling 
constant g used above T c , since we can also write k 2 = g/2f. The RG functions have no 
longer a polynomial form in k, being polynomial only in the gauge coupling constant /. The 
effective gauge coupling /(/«) determined by the condition j3f = shows that the charge is 
small for k sufficiently close to k = 1/ a/2 from above. Interestingly, the fixed point lies 
close to k = 1/ y/2. 

We have compared our values for the critical exponents with those found in the literature. 
Our one-loop value of u is larger than the XY value obtained from Monte Carlo simulations 
i"5fl and from RG analysis of the disorder field theory |12] . Using the expansion in powers of 



An of Section IV we have obtained r\ ~ —0.03 and v pa 0.92. With these values of rj and v 
we obtain a critical exponent 7 = v{2 — rj) ~ 1.87. This agrees with the Pade approximant 



analysis of the two-loop e-expansion obtained by Folk and Holovatch [19|. On the basis 
of the Pade approximants, these authors obtained the value v pa 0.857, showing a ~ 7% 
numerical difference from our result v ~ 0.92. The exponent 7 obtained from recent Monte 
Carlo simulations is 7 ~ 1.45 |33] . 



One question that immediately arises is how the critical behavior obtained with our 
method is related to approaches where T c is approached from above. In the pure O(N) 
symmetric 4 theory the critical singularities above and below T c are known to be the same. 
This result is obtained by analysing the Ward identities of the theory. For a dimensionally 
regularized theory, the same counterterms used in the disordered phase also renormalize 
the theory in the ordered phase fl|]. Since the critical exponents in the e-expansion are a 
direct consequence of the 1/e singularities included in the counterterms, we conclude that 
the critical exponents obtained by approaching T c from below are the same as when T c is 
approached from above. This should be true even if we work in fixed dimension. If we use 
the e-expansion in the GL model, the same argument holds. The problem in this case is that 
the superconducting fixed point appears only for sufficiently large N. We have seen that 
the value of N can be dramatically reduced if we consider a fixed dimension approach for 
T > T c . Unfortunately this is not enough to obtain a charged fixed point for N = 2. Then 
we have shown that by approaching T c from below we reach a charged fixed point for the 
physically interesting N = 2 case. It is not obvious which renormalization scheme above T c 
corresponds to our scheme below T c . Clearly, by working above T c we would never find the 
same renormalization constants. The point is that k, which appears explicitly below T c in a 
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non-polynomial form, can be only defined above T c through the ratio between the coupling 
constants. The vector potential is massive only below T c . 

We have discussed in detail the agreement of our theoretical ^-exponent with recent 
Monte Carlo simulations. It is now definite that it is negative. The only possible way to 
enforce a positive 77 > in the GL model is by adding an appropriate mass term, like the 



Proca term considered in Ref. EH. This should not be confused with our calculation of 
Section II where M — > at the end. A negative rj is caused by the attraction of the vortex 
lines on the average which causes a fluctuating vortex globule to have a smaller size than a 
free random chain. This point of view is related to the geometric interpretation considered 
in Ref. [H]. The addition of a Chern-Simons term which changes the size of vortex globules 
by a topological interaction seems incapable of producing a positive rj, since entanglement of 
vortex lines also tends to contract the globules. Indeed, one-loop calculation gives a negative 



77 for all values of the topological coupling [44, 45 



There is, however, the possibility that for larger topological coupling rj could become 
positive, since for infinite topological coupling 77 — > 0. This is suggested when calculating rj 
in a 1/N expansion [ d2[ . 

We hope that the new approach introduced in this paper will stimulate further discussions 
of the nature of the superconducting phase transition. In particular, we call for two- and 
higher- loop calculations to substantiate our claim that Kl — «*, and that the type-I-type-II 
crossover point is a Lifshitz point. 
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APPENDIX A: USEFUL INTEGRALS 

In this appendix we shall write some basic momentum space integrals in three dimensions 
that are used to obtain the results in the text. 
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The simplest integral, arising when computing tadpoles, is 



d 3 k 1 A m ,. s ,. , 

— + 0(m/A), (Al) 



(2tt) 3 k 2 + m 2 2tt 2 4tt 
where the ultraviolet cutoff is assumed to be very large, A>m. 

The following type of integral appears in the calculation of the two-point functions: 

/d 3 k 1 1 
7^~~vs 77 T72^ 2ui 2 i 2T = ^^ a (P' m i' m2 )' ( A2 ) 
(27r) 3 [(p — k) 2 + mfj(k 2 + ml) 87r|p| 

where a(p,m 1 ,m 2 ) is defined in Eq. ( |48|) . The following particular cases of the above 
integral are relevant: 

f d 3 k 1 1 

J (2tt) 3 [(p - k) 2 + m 2 ](k 2 + m 2 ) ~ (P ' m) ' ( } 

where 6(p, m) is defined in Eq. fl49f) , and 

^ 1 1 (A4) 



(27r) 3 (k 2 + mf ) (k 2 + mf) 47r(mi + m 2 ) ' 
When computing loop integrals involving the vector potential propagator we need the 

integrals: 

/ (2^ [( P -k) 2 + t 2 ](k 2 + m 2 ) = Th^ [2|pI (m2 - mi) + (p2 + ml ~ mlHp ' mi ' m2)] ' 

(A5) 
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(2vr) 3 [(p - k) 2 + m 2 ] (k 2 + m 2 ) 32vr 



mi + 777,2 — 



• [P +K+^)J[P +^-^)\ (p , mi , m2 )^ + 7 g^ {3m2 -5mi 



2P 



32tt P 2 



3(mi + m 2 )(mi — m 2 ) 2 3(p 2 + m?) 2 — 6m?m| + 3m| — 2m|p 2 1 
^ 1 ^173 a(p,mi,m 2 ) |(A6) 
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XX 



FIG. 1: Feynman graph contributing to the four-point function. The external lines correspond to 
the scalar fields, while the wiggles are the vector potential propagators. 




FIG. 2: Graphs contributing to the evaluation of the ir (a) and vector potential (b) two-point 
functions. The solid lines indicate the cr-propagators, the dashed ones 7r-propagators. The wiggles 
represent the vector potential. 
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FIG. 3: Behavior of effective square charge /(«). To the left of the vertical dashed line marking 
the the interval < y/2n < 0.096, J(k) is very large and positive, starting out from /(0) « 53.3 
and going to infinity for k — > 0.096/\/2. 
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